home *** CD-ROM | disk | FTP | other *** search
/ Mac Format 1995 June / MacFormat 25.iso / Shareware City / Developers / fortran-to-c-translator-11 / Mac F2C 1.1 / Test Project ƒ / test.cp (C++ Output) < prev    next >
Text File  |  1995-01-13  |  16KB  |  532 lines
  1. /* test.f -- translated by f2c (version 19941113).
  2.    You must link the resulting object file with the libraries:
  3.     -lf2c -lm   (in that order)
  4. */
  5.  
  6. #ifdef __cplusplus
  7. extern "C" {
  8. #endif
  9. #include "f2c.h"
  10.  
  11. /* Table of constant values */
  12.  
  13. static integer c__9 = 9;
  14. static integer c__1 = 1;
  15. static integer c__10 = 10;
  16. static integer c__3 = 3;
  17. static integer c__4 = 4;
  18. static integer c__5 = 5;
  19.  
  20.  
  21. /* Main program */ MAIN__()
  22. {
  23.     /* Format strings */
  24.     static char fmt_99[] = "(a1)";
  25.  
  26.     /* Builtin functions */
  27.     integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
  28.         e_wsle(), s_rsfe(cilist *), do_fio(integer *, char *, ftnlen), 
  29.         e_rsfe();
  30.     /* Subroutine */ int s_stop(char *, ftnlen);
  31.  
  32.     /* Local variables */
  33.     char junk[2];
  34.     extern /* Subroutine */ int i_o_test__(), flt_test__(integer *), 
  35.         int_test__(integer *), trn_test__();
  36.  
  37.     /* Fortran I/O blocks */
  38.     static cilist io___1 = { 0, 6, 0, 0, 0 };
  39.     static cilist io___2 = { 0, 6, 0, 0, 0 };
  40.     static cilist io___3 = { 0, 5, 0, fmt_99, 0 };
  41.     static cilist io___5 = { 0, 6, 0, 0, 0 };
  42.     static cilist io___6 = { 0, 6, 0, 0, 0 };
  43.     static cilist io___7 = { 0, 5, 0, fmt_99, 0 };
  44.     static cilist io___8 = { 0, 6, 0, 0, 0 };
  45.     static cilist io___9 = { 0, 6, 0, 0, 0 };
  46.     static cilist io___10 = { 0, 5, 0, fmt_99, 0 };
  47.     static cilist io___11 = { 0, 6, 0, 0, 0 };
  48.     static cilist io___12 = { 0, 6, 0, 0, 0 };
  49.     static cilist io___13 = { 0, 5, 0, fmt_99, 0 };
  50.     static cilist io___14 = { 0, 6, 0, 0, 0 };
  51.     static cilist io___15 = { 0, 6, 0, 0, 0 };
  52.     static cilist io___16 = { 0, 6, 0, 0, 0 };
  53.     static cilist io___17 = { 0, 6, 0, 0, 0 };
  54.     static cilist io___18 = { 0, 6, 0, 0, 0 };
  55.  
  56.  
  57. /*     This is a FORTRAN program to test Mac F2C v1.1 */
  58.     s_wsle(&io___1);
  59.     do_lio(&c__9, &c__1, "*****   Input/Output Test   *****", 33L);
  60.     e_wsle();
  61.     i_o_test__();
  62.     s_wsle(&io___2);
  63.     do_lio(&c__9, &c__1, "\n*****   End of I/O test, hit return to continue.\
  64. ..", 51L);
  65.     e_wsle();
  66.     s_rsfe(&io___3);
  67.     do_fio(&c__1, junk, 2L);
  68.     e_rsfe();
  69.     s_wsle(&io___5);
  70.     do_lio(&c__9, &c__1, "\n*****   Integer Math Test   *****", 34L);
  71.     e_wsle();
  72.     int_test__(&c__10);
  73.     s_wsle(&io___6);
  74.     do_lio(&c__9, &c__1, "\n*****   End of integer math test, hit return to \
  75. continue...", 60L);
  76.     e_wsle();
  77.     s_rsfe(&io___7);
  78.     do_fio(&c__1, junk, 2L);
  79.     e_rsfe();
  80.     s_wsle(&io___8);
  81.     do_lio(&c__9, &c__1, "\n*****   Floating Point Math Test   *****", 41L);
  82.     e_wsle();
  83.     flt_test__(&c__10);
  84.     s_wsle(&io___9);
  85.     do_lio(&c__9, &c__1, "\n*****   End of floating point math test, hit ret\
  86. urn to continue...", 67L);
  87.     e_wsle();
  88.     s_rsfe(&io___10);
  89.     do_fio(&c__1, junk, 2L);
  90.     e_rsfe();
  91.     s_wsle(&io___11);
  92.     do_lio(&c__9, &c__1, "\n*****   Transcendental Function Test   *****", 
  93.         45L);
  94.     e_wsle();
  95.     trn_test__();
  96.     s_wsle(&io___12);
  97.     do_lio(&c__9, &c__1, "\n*****   End of transcendental function test, hit\
  98.  return to continue...", 71L);
  99.     e_wsle();
  100.     s_rsfe(&io___13);
  101.     do_fio(&c__1, junk, 2L);
  102.     e_rsfe();
  103.     s_wsle(&io___14);
  104.     do_lio(&c__9, &c__1, "##################################################\
  105. ########################", 74L);
  106.     e_wsle();
  107.     s_wsle(&io___15);
  108.     do_lio(&c__9, &c__1, " If you noticed that floating point values did not\
  109.  round correctly when", 71L);
  110.     e_wsle();
  111.     s_wsle(&io___16);
  112.     do_lio(&c__9, &c__1, " displayed, please read the enclosed file \"If Flo\
  113. ats Don't Display Right\"", 73L);
  114.     e_wsle();
  115.     s_wsle(&io___17);
  116.     do_lio(&c__9, &c__1, "##################################################\
  117. ########################", 74L);
  118.     e_wsle();
  119.     s_wsle(&io___18);
  120.     do_lio(&c__9, &c__1, "\n*****   This completes all of the tests   *****", 
  121.         48L);
  122.     e_wsle();
  123.     s_stop("", 0L);
  124.     return 0;
  125. } /* MAIN__ */
  126.  
  127. /**************************************************************************/
  128.  
  129. /*     Subroutine to do the I/O tests */
  130.  
  131. /**************************************************************************/
  132. /* Subroutine */ int i_o_test__()
  133. {
  134.     /* Format strings */
  135.     static char fmt_399[] = "(a1)";
  136.     static char fmt_304[] = "(5x,a20,5(i1,2x),5x,5(f4.2,2x))";
  137.  
  138.     /* System generated locals */
  139.     olist o__1;
  140.     cllist cl__1;
  141.  
  142.     /* Builtin functions */
  143.     integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
  144.         e_wsle(), s_rsle(cilist *), e_rsle(), s_rsfe(cilist *), do_fio(
  145.         integer *, char *, ftnlen), e_rsfe();
  146.     void s_copy(char *, char *, ftnlen, ftnlen);
  147.     integer f_open(olist *), s_wsue(cilist *), do_uio(integer *, char *, 
  148.         ftnlen), e_wsue(), f_clos(cllist *), s_wsfe(cilist *), e_wsfe(), 
  149.         s_rsue(cilist *), e_rsue();
  150.  
  151.     /* Local variables */
  152.     char text[40];
  153.     real a[5];
  154.     integer i, j[5];
  155.     real x;
  156.     doublereal dx;
  157.  
  158.     /* Fortran I/O blocks */
  159.     static cilist io___19 = { 0, 6, 0, 0, 0 };
  160.     static cilist io___20 = { 0, 5, 0, 0, 0 };
  161.     static cilist io___22 = { 0, 6, 0, 0, 0 };
  162.     static cilist io___23 = { 0, 6, 0, 0, 0 };
  163.     static cilist io___24 = { 0, 5, 0, 0, 0 };
  164.     static cilist io___26 = { 0, 6, 0, 0, 0 };
  165.     static cilist io___27 = { 0, 6, 0, 0, 0 };
  166.     static cilist io___28 = { 0, 5, 0, 0, 0 };
  167.     static cilist io___30 = { 0, 6, 0, 0, 0 };
  168.     static cilist io___31 = { 0, 6, 0, 0, 0 };
  169.     static cilist io___32 = { 0, 5, 0, 0, 0 };
  170.     static cilist io___34 = { 0, 6, 0, 0, 0 };
  171.     static cilist io___35 = { 0, 6, 0, 0, 0 };
  172.     static cilist io___36 = { 0, 5, 0, fmt_399, 0 };
  173.     static cilist io___39 = { 0, 60, 0, 0, 0 };
  174.     static cilist io___40 = { 0, 6, 0, 0, 0 };
  175.     static cilist io___41 = { 0, 6, 0, fmt_304, 0 };
  176.     static cilist io___42 = { 0, 50, 0, 0, 0 };
  177.     static cilist io___43 = { 0, 6, 0, 0, 0 };
  178.     static cilist io___44 = { 0, 6, 0, fmt_304, 0 };
  179.  
  180.  
  181. /*     Screen I/O tests */
  182.     s_wsle(&io___19);
  183.     do_lio(&c__9, &c__1, "\nPart 1:  Screen I/O tests.\n\nEnter an integer v\
  184. alue.", 52L);
  185.     e_wsle();
  186.     s_rsle(&io___20);
  187.     do_lio(&c__3, &c__1, (char *)&i, (ftnlen)sizeof(integer));
  188.     e_rsle();
  189.     s_wsle(&io___22);
  190.     do_lio(&c__9, &c__1, "The number you entered was:", 27L);
  191.     do_lio(&c__3, &c__1, (char *)&i, (ftnlen)sizeof(integer));
  192.     e_wsle();
  193.     s_wsle(&io___23);
  194.     do_lio(&c__9, &c__1, "\nEnter a single precision floating point value...",
  195.          49L);
  196.     e_wsle();
  197.     s_rsle(&io___24);
  198.     do_lio(&c__4, &c__1, (char *)&x, (ftnlen)sizeof(real));
  199.     e_rsle();
  200.     s_wsle(&io___26);
  201.     do_lio(&c__9, &c__1, "The number you entered was: ", 28L);
  202.     do_lio(&c__4, &c__1, (char *)&x, (ftnlen)sizeof(real));
  203.     e_wsle();
  204.     s_wsle(&io___27);
  205.     do_lio(&c__9, &c__1, "\nEnter a double precision floating point value...",
  206.          49L);
  207.     e_wsle();
  208.     s_rsle(&io___28);
  209.     do_lio(&c__5, &c__1, (char *)&dx, (ftnlen)sizeof(doublereal));
  210.     e_rsle();
  211.     s_wsle(&io___30);
  212.     do_lio(&c__9, &c__1, "The number you entered was: ", 28L);
  213.     do_lio(&c__5, &c__1, (char *)&dx, (ftnlen)sizeof(doublereal));
  214.     e_wsle();
  215.     s_wsle(&io___31);
  216.     do_lio(&c__9, &c__1, "\nEnter some text (40 char max)...", 33L);
  217.     e_wsle();
  218.     s_rsle(&io___32);
  219.     do_lio(&c__9, &c__1, text, 40L);
  220.     e_rsle();
  221.     s_wsle(&io___34);
  222.     do_lio(&c__9, &c__1, "The text you entered was: ", 26L);
  223.     do_lio(&c__9, &c__1, text, 40L);
  224.     e_wsle();
  225.     s_wsle(&io___35);
  226.     do_lio(&c__9, &c__1, "\nPart 2:  file I/O tests.  Hit return to continue\
  227. ...", 52L);
  228.     e_wsle();
  229.     s_rsfe(&io___36);
  230.     do_fio(&c__1, text, 40L);
  231.     e_rsfe();
  232. /*     File I/O tests:  Store some values and write them to file */
  233.     for (i = 1; i <= 5; ++i) {
  234.     j[i - 1] = i;
  235.     a[i - 1] = (doublereal) i;
  236.     }
  237.     s_copy(text, "A test message.", 40L, 15L);
  238.     o__1.oerr = 0;
  239.     o__1.ounit = 60;
  240.     o__1.ofnmlen = 8;
  241.     o__1.ofnm = "test.dat";
  242.     o__1.orl = 0;
  243.     o__1.osta = 0;
  244.     o__1.oacc = 0;
  245.     o__1.ofm = "unformatted";
  246.     o__1.oblnk = 0;
  247.     f_open(&o__1);
  248.     s_wsue(&io___39);
  249.     do_uio(&c__1, text, 40L);
  250.     do_uio(&c__5, (char *)&j[0], (ftnlen)sizeof(integer));
  251.     do_uio(&c__5, (char *)&a[0], (ftnlen)sizeof(real));
  252.     e_wsue();
  253.     cl__1.cerr = 0;
  254.     cl__1.cunit = 60;
  255.     cl__1.csta = 0;
  256.     f_clos(&cl__1);
  257.     s_wsle(&io___40);
  258.     do_lio(&c__9, &c__1, "Wrote the following data to file test.dat:\n", 43L);
  259.     e_wsle();
  260.     s_wsfe(&io___41);
  261.     do_fio(&c__1, text, 40L);
  262.     for (i = 1; i <= 5; ++i) {
  263.     do_fio(&c__1, (char *)&j[i - 1], (ftnlen)sizeof(integer));
  264.     }
  265.     for (i = 1; i <= 5; ++i) {
  266.     do_fio(&c__1, (char *)&a[i - 1], (ftnlen)sizeof(real));
  267.     }
  268.     e_wsfe();
  269. /* Reset the variables and read them back */
  270.     for (i = 1; i <= 5; ++i) {
  271.     j[i - 1] = 99;
  272.     a[i - 1] = (float)99.;
  273.     }
  274.     s_copy(text, "reset", 40L, 5L);
  275.     o__1.oerr = 0;
  276.     o__1.ounit = 50;
  277.     o__1.ofnmlen = 8;
  278.     o__1.ofnm = "test.dat";
  279.     o__1.orl = 0;
  280.     o__1.osta = 0;
  281.     o__1.oacc = 0;
  282.     o__1.ofm = "unformatted";
  283.     o__1.oblnk = 0;
  284.     f_open(&o__1);
  285.     s_rsue(&io___42);
  286.     do_uio(&c__1, text, 40L);
  287.     do_uio(&c__5, (char *)&j[0], (ftnlen)sizeof(integer));
  288.     do_uio(&c__5, (char *)&a[0], (ftnlen)sizeof(real));
  289.     e_rsue();
  290.     cl__1.cerr = 0;
  291.     cl__1.cunit = 50;
  292.     cl__1.csta = 0;
  293.     f_clos(&cl__1);
  294.     s_wsle(&io___43);
  295.     do_lio(&c__9, &c__1, "\nRead the following data from file test.dat:\n", 
  296.         45L);
  297.     e_wsle();
  298.     s_wsfe(&io___44);
  299.     do_fio(&c__1, text, 40L);
  300.     for (i = 1; i <= 5; ++i) {
  301.     do_fio(&c__1, (char *)&j[i - 1], (ftnlen)sizeof(integer));
  302.     }
  303.     for (i = 1; i <= 5; ++i) {
  304.     do_fio(&c__1, (char *)&a[i - 1], (ftnlen)sizeof(real));
  305.     }
  306.     e_wsfe();
  307.     return 0;
  308. } /* i_o_test__ */
  309.  
  310. /**************************************************************************/
  311.  
  312. /*     Subroutine to do the integer math tests */
  313.  
  314. /**************************************************************************/
  315. /* Subroutine */ int int_test__(integer *m)
  316. {
  317.     /* Format strings */
  318.     static char fmt_203[] = "(10x,\002n\002,5x,\002n^2\002,5x,\002n^3\002,\
  319. 5x,\002n/2\002,3x,\002n^2/2\002,3x,\002n^3/2\002)";
  320.     static char fmt_202[] = "(5x,6(i6,2x))";
  321.  
  322.     /* System generated locals */
  323.     integer i__1, i__2, i__3, i__4;
  324.  
  325.     /* Builtin functions */
  326.     integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
  327.         e_wsle(), s_wsfe(cilist *), e_wsfe(), do_fio(integer *, char *, 
  328.         ftnlen);
  329.  
  330.     /* Local variables */
  331.     integer i, j, k;
  332.  
  333.     /* Fortran I/O blocks */
  334.     static cilist io___45 = { 0, 6, 0, 0, 0 };
  335.     static cilist io___46 = { 0, 6, 0, fmt_203, 0 };
  336.     static cilist io___50 = { 0, 6, 0, fmt_202, 0 };
  337.  
  338.  
  339.     s_wsle(&io___45);
  340.     do_lio(&c__9, &c__1, "\nGenerate a table of integers, squares, cubes, an\
  341. d their halves.\n", 65L);
  342.     e_wsle();
  343.     s_wsfe(&io___46);
  344.     e_wsfe();
  345.     i__1 = *m;
  346.     for (i = 1; i <= i__1; ++i) {
  347. /* Computing 2nd power */
  348.     i__2 = i;
  349.     j = i__2 * i__2;
  350. /* Computing 3rd power */
  351.     i__2 = i, i__3 = i__2;
  352.     k = i__3 * (i__2 * i__2);
  353.     s_wsfe(&io___50);
  354.     do_fio(&c__1, (char *)&i, (ftnlen)sizeof(integer));
  355.     do_fio(&c__1, (char *)&j, (ftnlen)sizeof(integer));
  356.     do_fio(&c__1, (char *)&k, (ftnlen)sizeof(integer));
  357.     i__2 = i / 2;
  358.     do_fio(&c__1, (char *)&i__2, (ftnlen)sizeof(integer));
  359.     i__3 = j / 2;
  360.     do_fio(&c__1, (char *)&i__3, (ftnlen)sizeof(integer));
  361.     i__4 = k / 2;
  362.     do_fio(&c__1, (char *)&i__4, (ftnlen)sizeof(integer));
  363.     e_wsfe();
  364.     }
  365.     return 0;
  366. } /* int_test__ */
  367.  
  368. /**************************************************************************/
  369.  
  370. /*     Subroutine to do the floating point math tests */
  371.  
  372. /**************************************************************************/
  373. /* Subroutine */ int flt_test__(integer *m)
  374. {
  375.     /* Format strings */
  376.     static char fmt_205[] = "(12x,\002x\002,6x,\002x^2\002,6x,\002x^3\002,\
  377. 6x,\002x/2\002,4x,\002x^2/2\002,4x,\002x^3/2\002)";
  378.     static char fmt_201[] = "(5x,6(f8.2,1x))";
  379.  
  380.     /* System generated locals */
  381.     integer i__1;
  382.     real r__1, r__2, r__3;
  383.  
  384.     /* Builtin functions */
  385.     integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
  386.         e_wsle(), s_wsfe(cilist *), e_wsfe(), do_fio(integer *, char *, 
  387.         ftnlen);
  388.  
  389.     /* Local variables */
  390.     integer i;
  391.     real x1, x2, x3;
  392.  
  393.     /* Fortran I/O blocks */
  394.     static cilist io___51 = { 0, 6, 0, 0, 0 };
  395.     static cilist io___52 = { 0, 6, 0, fmt_205, 0 };
  396.     static cilist io___57 = { 0, 6, 0, fmt_201, 0 };
  397.  
  398.  
  399.     s_wsle(&io___51);
  400.     do_lio(&c__9, &c__1, "\nGenerate a table of floats, their squares, cubes\
  401. , and their halves.\n", 69L);
  402.     e_wsle();
  403.     s_wsfe(&io___52);
  404.     e_wsfe();
  405.     i__1 = *m;
  406.     for (i = 1; i <= i__1; ++i) {
  407.     x1 = i * (float)1.;
  408. /* Computing 2nd power */
  409.     r__1 = x1;
  410.     x2 = r__1 * r__1;
  411. /* Computing 3rd power */
  412.     r__1 = x1, r__2 = r__1;
  413.     x3 = r__2 * (r__1 * r__1);
  414.     s_wsfe(&io___57);
  415.     do_fio(&c__1, (char *)&x1, (ftnlen)sizeof(real));
  416.     do_fio(&c__1, (char *)&x2, (ftnlen)sizeof(real));
  417.     do_fio(&c__1, (char *)&x3, (ftnlen)sizeof(real));
  418.     r__1 = x1 / 2;
  419.     do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
  420.     r__2 = x2 / 2;
  421.     do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
  422.     r__3 = x3 / 2;
  423.     do_fio(&c__1, (char *)&r__3, (ftnlen)sizeof(real));
  424.     e_wsfe();
  425.     }
  426.     return 0;
  427. } /* flt_test__ */
  428.  
  429. /**************************************************************************/
  430.  
  431. /*     Subroutine to do the transcendental function tests */
  432.  
  433. /**************************************************************************/
  434. /* Subroutine */ int trn_test__()
  435. {
  436.     /* Format strings */
  437.     static char fmt_207[] = "(9x,\002x\002,10x,\002sin(x)\002,8x,\002cos(x\
  438. )\002,3x,\002sin(x)^2 + cos(x)^2\002)";
  439.     static char fmt_200[] = "(5x,i2,\002*pi/6\0023x,f11.8,3x,f11.8,3x,f15.10)"
  440.         ;
  441.     static char fmt_299[] = "(a1)";
  442.     static char fmt_208[] = "(11x,\002x\002,16x,\002log(x)\002,9x,\002exp(lo\
  443. g(x))\002)";
  444.     static char fmt_201[] = "(5x,f13.10,5x,f13.10,5x,f13.10)";
  445.  
  446.     /* System generated locals */
  447.     doublereal d__1;
  448.  
  449.     /* Builtin functions */
  450.     integer s_wsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
  451.         e_wsle(), s_wsfe(cilist *), e_wsfe();
  452.     double sin(doublereal), cos(doublereal);
  453.     integer do_fio(integer *, char *, ftnlen), s_rsfe(cilist *), e_rsfe();
  454.     double log(doublereal), exp(doublereal);
  455.  
  456.     /* Local variables */
  457.     char junk[2];
  458.     doublereal c;
  459.     integer i;
  460.     doublereal s, x, c2, s2, pi;
  461.  
  462.     /* Fortran I/O blocks */
  463.     static cilist io___59 = { 0, 6, 0, 0, 0 };
  464.     static cilist io___60 = { 0, 6, 0, 0, 0 };
  465.     static cilist io___61 = { 0, 6, 0, fmt_207, 0 };
  466.     static cilist io___68 = { 0, 6, 0, fmt_200, 0 };
  467.     static cilist io___69 = { 0, 6, 0, 0, 0 };
  468.     static cilist io___70 = { 0, 5, 0, fmt_299, 0 };
  469.     static cilist io___72 = { 0, 6, 0, 0, 0 };
  470.     static cilist io___73 = { 0, 6, 0, fmt_208, 0 };
  471.     static cilist io___74 = { 0, 6, 0, fmt_201, 0 };
  472.  
  473.  
  474.     pi = (float)3.141592653589793;
  475.     s_wsle(&io___59);
  476.     do_lio(&c__9, &c__1, "\nPart 1: Trig Functions", 23L);
  477.     e_wsle();
  478.     s_wsle(&io___60);
  479.     do_lio(&c__9, &c__1, "\nGenerate a table of x, sin(x), cos(x) and the su\
  480. m of their squares.\n", 69L);
  481.     e_wsle();
  482.     s_wsfe(&io___61);
  483.     e_wsfe();
  484.     for (i = 1; i <= 12; ++i) {
  485.     x = i * pi / (float)6.;
  486.     s = sin(x);
  487.     c = cos(x);
  488. /* Computing 2nd power */
  489.     d__1 = s;
  490.     s2 = d__1 * d__1;
  491. /* Computing 2nd power */
  492.     d__1 = c;
  493.     c2 = d__1 * d__1;
  494.     s_wsfe(&io___68);
  495.     do_fio(&c__1, (char *)&i, (ftnlen)sizeof(integer));
  496.     do_fio(&c__1, (char *)&s, (ftnlen)sizeof(doublereal));
  497.     do_fio(&c__1, (char *)&c, (ftnlen)sizeof(doublereal));
  498.     d__1 = s2 + c2;
  499.     do_fio(&c__1, (char *)&d__1, (ftnlen)sizeof(doublereal));
  500.     e_wsfe();
  501.     }
  502.     s_wsle(&io___69);
  503.     do_lio(&c__9, &c__1, "\nPart 2:  Exponential functions; hit return to co\
  504. ntinue...", 58L);
  505.     e_wsle();
  506.     s_rsfe(&io___70);
  507.     do_fio(&c__1, junk, 2L);
  508.     e_rsfe();
  509.     s_wsle(&io___72);
  510.     do_lio(&c__9, &c__1, "Generate a table of x, log(x), and exp(log(x))\n", 
  511.         47L);
  512.     e_wsle();
  513.     s_wsfe(&io___73);
  514.     e_wsfe();
  515.     for (i = 1; i <= 10; ++i) {
  516.     x = (doublereal) i;
  517.     s = log(x);
  518.     c = exp(s);
  519.     s_wsfe(&io___74);
  520.     do_fio(&c__1, (char *)&x, (ftnlen)sizeof(doublereal));
  521.     do_fio(&c__1, (char *)&s, (ftnlen)sizeof(doublereal));
  522.     do_fio(&c__1, (char *)&c, (ftnlen)sizeof(doublereal));
  523.     e_wsfe();
  524.     }
  525.     return 0;
  526. } /* trn_test__ */
  527.  
  528. /* Main program alias */ int test_f2c__ () { MAIN__ (); return 0; }
  529. #ifdef __cplusplus
  530.     }
  531. #endif
  532.